Term reference

Term Definition
Rs_annual annual soil respiration (Rs, g C m-2 yr-1)
Rs_annual_bahn annual Rs computed using Bahn (2010) method (g C m-2 yr-1)
amat annual mean air temperature (e.g., amat of 2000 is the evarage of 12 months’ air temperature in 2000)
mat mean annual air temperature within a long period (in this study is from 1964 to 2017)
amst annual mean soil temperature

Summary

Introduction

Rs measured at mean annual soil temperature

Rs measured at mean annual soil temperature

The objects of this analysis are

3. Methods

Data

Statistics

Update Bahn model

Update Bahn model: * new_model1: only update parameters, but same model formulation * new_model2: add other parameters to the model for better prediction

[JJTODO: consider whether we want to use “amst” or “mast”. Give the prevalance of “MAT” I feel like the latter would be clearer?]

We examined many possibilities for why the Rs_annual_bahn vs Rs_annual relationship is not 1:1.

3.4.1 Ts sources (MGRsD, MGRsD_TAIR, From paper, Rs_Ts_relationship)

First, we tested the effect of different soil temperature sources on the Rs_annual_bahn vs Rs_annual relationship:

  • From MGRsD means mean annual soil temperature (amst) are from a global monthly soil respiration database, each site has more than 12 months measured soil temperature read from original papers.
  • From paper means amst were reported from the original paper (table, figures, or description).
  • Partly from TAIR means: some studies did not measure soil temperature all year, for those months, we predict soil temperature based on monthly air temperature (Tsoil = 2.918 + 0.829*Tair). This model was developed based on the sites which have >= 12 months soil temperature measurements.

There are 67 records for which I cannot get the soil temperature information through the above three methods. In these cases, based on the Rs_Ts_relationship and reported Rs_annual, I calculated the amst.

  • The calculated amst was then compared with the annual Tair, if they are well matched (error < 5%), calculated mast were used.
  • When calculated amst and annual Tair do not match, this usually indicates a potential problem, and then I went back to the manuscript and checked.
  • Whenever a paper reported annual mean Ts, I compared the reported mast and estimated mast based on the Rs_Ts_relationship, and found they are well matched.

Generally, Ts sources do not have clear effects on the Rs_annual_bahn and Rs_annual relationship.

3.4.2 Annual Rs or Ts coverage effect

Second, we tested whether Ts and Rs coverage (e.g., 0-0.5 means Rs or Ts only measured less than 6 months, versus the entire year), but found that Ts and Rs coverage do not have significant effects on the Rs_annual_bahn vs Rs_annual relationship.

3.4.3 Effect of maximum allowed divergence between global climate data set and site-specific air temperature

Third, we tested whether maximum allowed divergence between global climate data set and site-specific air temperature affect the Rs_annual_bahn vs Rs_annual relationship.

As we throw out data points with high divergence, R2 and RSE go up and down inconsistently and by small amounts, suggesting that the Tair divergence do not have a consistent effect.

3.4.4 Effect of maximum allowed divergence between annual precipitation from paper and Del

Fourth, we also tested whether the maximum allowed divergence between the global climate data set and site-specific precipitation affects the Rs_annual_bahn vs Rs_annual relationship. In other words, does a bias in the global data affect things?

As we throw out data points with high divergence, R2 and RSE showed no large or consistent changes, suggesting that the precipitation divergence does not have a large effect.

We also compared: * MAT from U. Delaware (MAT_Del) and MAT reported from the papers (a) * TAnnual from U. Delaware (TAnnual_Del) and annual temperature from papers (study_temp) (b) * MAP from U. Delaware (MAP_Del) and MAP reported from the papers (c) * PAnnual from U. Delaware (PAnnual_Del) and annual precipitation from papers (study_precip) (d)

In general, the temperature and precipitation from the University of Delaware climate data matched the data reported from publications well. This supports the idea that the divergence between global climate data set and site-specific precipitation/temperature has little to no effect.

We tested the effect of precipitation and temperature variability (quantified by standard deviation from 1961 to 2014), using multiple linear regression, with divergence as catergorical indicator. [[JJTODO: this needs to be clearer. sd of what exactly? Annual means and sums?]]

We found that they have no significant effect on the Rs_annual_bahn vs Rs_annual relationship.

3.4.5 Test ecosystem type

The Rs_annual_ban vs Rs_annual relationship varies among different ecosystems.

  • For example, agriculture has lower slope but wetland has higher slope.
  • However, it is unlikely that the data from agriculture and wetland shift the overall regression between Rs_annual_bahn vs Rs_annual away from 1:1 line.
  • When we remove the Ag data, the Rs_annual_bahn vs Rs_annual regression still differs from 1:1 line.

  • When we plot Ag alone, we see the Rs_annual_bahn vs Rs_annual regression in Ag does not greatly differ from the rest.
  • We also examined the influential outlier points (cooks.distance > 0.5). When the outliers were removed, the regression showed no difference, indicating that outliers do not have large effects.
  • Similar conclusion for wetland data.

3.4.6 Test Rs measurement method

Different measurement methods do not affect the Rs_annual_bahn vs Rs_annual relationship.

3.4.7 RA- or RH-dominated effect?

A particularly interesting question is whether the Rs_annual_bahn vs Rs_annual relationship changes in sites dominated by autotrophic (RA) or heterotrophic (RH) respiration. This might be the case if, for example, one respiration source had a consistently highly temperature sensitivity.

  • RA dominated sites tend to have larger intercept than RH dominated sites, but no difference in slope.

We tested Q10 (temperature sensitivity) and R10 (Rs at a standardized 10 C) at RA and RH-dominated sites.

  • RA dominated sites have larger Q10 (over the 0-10, 5-10, 10-20, and 0-20 soil temperature ranges) values and R10 values.

3.4.8 Biome effect?

  • ‘Mediterranean’ sites exhibit large differences from other biomes.

3.4.9 Drought effect?

In their paper, Bahn et al. (2010) reported that drought stress significantly affected the relationship. Using our new datasets we found that:

  • MAP effects on the Rs_annual_bahn vs Rs_annual relationship is very limited (p ~~ 0.05).
  • The slope and intercept changes do not followed a clear pattern as MAP changes.
  • MAP is not a good drought index.

We then tested standardized drought index (SPI).

  • SPI significantly affects Rs_annual_bahn vs Rs_annual relationship; as SPI decreases, the slope decreases, means the Bahn approach tends to overestimate Rs_annual under drought condition.

Since SPI compares annual precipitation at a site with average precipitation over a period (we used 1964-2014), it can not describe spatial drought in drought.

We thus also used another drought index, the Palmer Drought Severity Index (PDSI), to characterize spatial effects.

  • PDSI significantly affects the Rs_annual_bahn vs Rs_annual relationship: as it becomes drier, the slope tend to decrease, meaning that the Bahn method overpredicts the annual flux. This is consistent with the SPI finding above.

4 Update Bhan model

The previous analysis show that Rs_annual_bahn does not well represent Rs-annual; we tested several possibilities to understand why, but no solutions were found to make the original model predict Rs_annual robustly.

It is possible that the Bahn (2010) model only used 80 sites across globel, it is not representive. * We thus updated the Bahn (2010) model’s parameters (but with same formulation, named new1 model). * Following Bahn (2010), and because of the test performed above, we built a model for Mediterranean (n=21), and another model for the rest of the data (n=802). [JJTODO: calculate]

[JJTODO: the graph below needs a legend–very unclear without one]

[JJTODO: report model stats/parameters?]

[JJTODO: need a summary statement here.]

4.2 Predicting Rs_annual from air temperatures

Since high resolution soil temperature is still lacking, and/or has lower accuracy than air temperature data, we want to test whether we can use Rs at annual mean air temperature (amat) or mean annual temperature (mat) to predict Rs_annual. [JJTODO: clarify the difference between these]

  • The regression between Rs_annual_bahn_amat and Rs_annual falls away from the 1:1 line, with an intercept significantly different from 0 (p<0.001) and slope significantly different from 1 (p<0.001).
  • Note that we tested both new1 and new2 models.
  • Using amat or mat show no big difference, but it is surprising that using mat is slightly better than using T_Annual. [JJTODO: be consistent in terminology]

We detected 2 outliers in the Rs_annual_bahn vs Rs_annual regression. * Remove these two outliers significantly improved the model (slope changed from 0.78 to 0.87, intercept decreased from 222 to 157, however, p values for slope and intercept are still < 0.001).

## Warning: Removed 2 rows containing non-finite values (stat_smooth).
## Warning: Removed 2 rows containing missing values (geom_point).

## Warning: Removed 1 rows containing non-finite values (stat_smooth).
## Warning: Removed 1 rows containing missing values (geom_point).

[JJTODO: this is very unclear] * We adjusted the air temperature by TAnnual_adj = 2.918+0.829TAnnual MAT_adj = 2.918+0.829MAT The results are better, but still not resolve the problem (p<0.05)

4.3 Re-simulate a model (new3)

We re-calculated soil respiration at annual mean air temperature (Rs_amat, i.e., using air temperature rather than soil temperature to calculate soil respiration).

  • Then we re-simulated the relationship between Rs_annual and Rs_amat.
  • Rs_annual = 729.09225 * (Rs_amat ^ 0.46535) + 89.7789 * spi – Medeterrean
  • Rs_annual = 588.618 * ( Rs_amat ^ 0.65022 ) + 22.59026 * spi + 11.29775*pdsi – exclude Mediterranean sites
  • If we update the model, Rs_annual_bahn_amat can represent Rs_annual.

5. Discussion & questions

These results have a direct bearing on two important problems for Rs and more generally carbon-cycle measurement and modeling: * We have many more measurements Rs in mid-latitude regions and developed countries.Less-developed countries are constrained by lack of resources, and thus we do not have enough measurements from spouth hetmesphere, arctic, and tropical regions (Xu and Shang 2016) * It is difficult to measure soil respiration all year around in cold regions, but critical because of high rates of climate change and large soil C stocks

Global spatial distribution of soil respiration sites

Global spatial distribution of soil respiration sites

We show that Rs measured at annual mean temperature (soil temperature or air temperature) can represent Rs_annual well, with well-quantified errors. This capability could be used to improve Rs measure frequency and greatly decrease cost, which becomes more important in the southern hemisphere and cold regions.

6. More analysis in the future

7. Select graphs for poster

## Warning: Removed 1 rows containing non-finite values (stat_smooth).
## Warning: Removed 1 rows containing missing values (geom_point).
## Warning: Removed 1 rows containing non-finite values (stat_smooth).
## Warning: Removed 1 rows containing missing values (geom_point).
## Warning: Removed 1 rows containing non-finite values (stat_smooth).
## Warning: Removed 1 rows containing missing values (geom_point).
## Warning: Removed 1 rows containing non-finite values (stat_smooth).
## Warning: Removed 1 rows containing missing values (geom_point).